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We study the problem of a boundary magnetic impurity coupled with the solvable t — J chains. 
0^ ' Our model provides a good start point to understand the Kondo problem in a Luttinger liquid as well 

as in a strongly correlated host. The Kondo coupling constant J can take arbitrary values, which 
allows us to study the ferromagnetic and antiferromagnetic Kondo problems simultaneously. It is 
O |. shown that both the Kondo coupling and the scalar potential effectively reconcile the spin dynamics 

at low temperatures. Generally, the impurity spin is split into two ghost spins via coupling effect. 
The residual entropy, which can be exactly derived for the present model, is strongly interaction- 
I , dependent. This hints the local spin configuration is very complicated and very different from that 

in the conventional Kondo problem. An unscreening phenomenon in the antiferromagnetic regime 
is found, which reveals the impurity potential plays an important role for the Kondo problem in a 
. strongly correlated host. The low temperature specific heat is calculated in the framework of local 

,1 ^ ' Landau-Luttinger liquid theory. 
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JL ■ I. INTRODUCTION 

a ■ 

With the development of the nanofabrication techniques for quantum wires and the prediction of edge stateS|-i|H the 
quantum Hall effect, the interest in one-dimensional (ID) quantum systems has been renewed in recent yearsou. In 
fact, much of the interest in ID quantum systems is due to Anderson's observation!] that the normal state properties 
y—i , of the quasi 2D high-T c superconductors are strikingly different from all known metals and can not be reconciled with 
^ ■ Landau's Fermi liquid theory but are more similar to properties of ID metals. In another hand, the impurity problem 
has been a current interest in the field of condensed matter physics. A well known example is the Kondo problem, 
which stimulated a strong challenge to traditional perturbation theory and provided a possible "laboratory" to search 
for the non-Fermi liquid behavior. The local perturbation problem to a ID Fermi system has been the subject of 
an intensive theoretical investigation in the recent years, for both its interesting anomalies with respect to that of a 
higher-dimensiaa^il system, and its relevance to a variety physical situations such as the-.transport behavior jpf the 
^\ , quantum wiresBu and the tunneling through a constriction in the quantum Hall regimea. Kane and Fishcrtl have 
^ ■ argued that a single impurity in a ID repulsive interacting system in fact corresponds to a chain disconnected at the 
impurity site at low energy scales. Their observation was also justified in the framework of the renormalization group 
analysis of boundary conformal field theoryLri, which show that the open boundary condition is the stable fixed point 
for repulsive interaction. An interesting problem is the ID electron system couples with a magnetic impurity or the 
Kondo problem in a ID metals. This problem was first considered by Lee and TonertJ, who found the crossover of 
the Kondo temperature from power law dependence on the Kondo coupling constant to an exponential one, when 
q , the electron correlation-goes from strong limit to weak limit. Subsequently, a poor man's scaling was performed by 
Furusaki and NagaosaEj, who addressed a conjecture thatpierromagnetic Kondo screening may occur in ID due to 
the special topology. The boundary conformal field theoryliil gave out a classification of critical behavior for the ID 
Kondo problem (without impurity potential). It turns there are only two possibilities, a standard low-temperature 
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thermodynamics or a non-Fermi liquid observed by Furusaki and Nagaosa. The non- Fermi liquid behavior is induced 
by the tunneling effect of conduction electrons through the impurity which depends only on the bulk properties but 
not on the detail of the impurityta. 

Despite such important progress, the problem of few impurities (potential, magnetic, especially both) embedded 
in a strongly correlated ID system is still not very well understood. We remark the study on the integrable models 
generally gives some useful information to understand a variety of physical situations in ID, if not fully. There are a few 
integrable models related to the impurity problem in ID quantum svstems: an impurity spin embedded in a spin-1/2 
Heisenberg jChain solved many years ago by Andrei and JohannessonL3; an integrable impurity in the supersymmetric 
t — J modeltj solved by Bediirftig et al.. The impurities in these models are exactly transparent due to the unphysical 
terms in the hamiltonians ancL^aatradict to the predictions of Kane and Fisheiu. This shortcoming was overcome in 
some other integrable modelsEj~LJ by introducing the open boundary condition at the impurity site. Unlike in 3D, 
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the scalar potentials! and the bond deformatiorj0 have significant effects to the low-temperature thermodynamics of 
the impurity. The anomaly residual entropy induced by the boundary effect was predicted in these models, which 
strongly suggests that the charge-spin cooperation in ID correlated electron systems may play an important role. 

In this paper, we consider a boundary magnetic impurity coupled with the integrable t — J chains. This is the 
detailed and extended description of our earlier letterlia. The Hamiltonian we shall study reads 

H = H + Hi, 

L-l L-l 

H " = - t J2 53( c i,<rCi+i,<7 + h.c.) + J2( JS i ■ S J+i + Vn i n i+i)> (1) 

3=1 3=1 

Hi = 2J& L • S + V t n L , 

where L is the length or site number of the chain; c\ a (cj, a ) are the creation (annihilation) operators of the conduction 

electrons and Sj = 1/2 ^ CT a , cj a r at(T 'Cj tCr i are the spin operators of the conduction electrons; rij lCr = c ] a c j,° IS the 
number operator; S is the impurity spin with strength S; t, J, V, Ji and Vi are all constants. Notice that the single 
occupation condition nj < 1 is assumed in (1). Without losing generality, we shall put t = 1 in the following text. 

The structure of the present paper is the following: In the subsequent section- we derive the integrable conditions 
for a variety of choices of the parameters J and V. The reflection equationEj constrains the integrability of the 
present model. That means the parameters Ji and Vi can not be independent but are parameterized by an unique 
c-number c to ensure the integrability of the model. In addition, the Bethe ansatz equations for the integrable cases 
will be given. In sect. Ill, we study the ground state properties and low-lying excitations. The formation of boundary 
bound states will be discussed in detail. Sect. IV is attributed to the derivation of the thermodynamic equations. The 
residual entropy induced by the impurity will be discussed. A different method, i.e., the local Landau-Luttinger liquid 
descriptio nESE-D is used to study the low-temperature thermodynamics of the impurity. Such a method can also be 
applied to non-integrable models and is thus general to the impurity problem in the ID quantum systems. Concluding 
remarks will be given in Sect.V. In appendix A, the general solution for the reflection matrix is given. The appendix 
B is attributed to the eigenvalue problem of the nested Bethe ansatz. 

II. INTEGRABLE CONDITIONS 

It is well known that Hq is exactly solvable for J = ±2, V = =Fl/2, ±3/200. By including the impurity, any 
electron impinging on the impurity will be completely reflected and suffer a reflection matrix Rj. The waves are 
therefore reflected at either end as 

e tk ]X fl i (jb i ) e -*^»-!M»y£ ) x ~ L. (2) 
Let us consider the two particle case. There are two ways from an initial state (fci, k%, |) to a final state (— k%, — fe, |): 
I-(k 1 ,k 2 ,\) -> (k 2 ,k 1 ,\) -> (k 2 ,-k 1 ,\) -> (-fci,fc 2 ,|) -> (-fa, -k 2 , |), 

II.(ki,k 2 , |) -> {fa, -k 2 , |) -> (-k 2 ,fa, |) -> (-fe, -fa, |) -> (-fa, -ha, I), 

where the symbol | denotes the open boundary. Since the physical process is unique, the following equation must hold 

S 12 (fa, k 2 )Ri(fa)S 12 (fa,-k 2 )R 2 (k 2 ) = R 2 (k 2 )S 12 (fa, -k^R^S^fa, k 2 ). (3) 

Above S\ 2 is the scattering matrix between the two electrons. This is just the reflection equationEl. For the multi- 
particle cases, p. long as the scattering matrix is factorizable or the two-body scattering matrix satisfies the Yang- 
Baxter relatione^ 

Si 2 (fa, k 2 )Si 3 (fa, k 3 )S 23 (k 2 , fc 3 ) = S 23 (k 2 , k 3 )Si 3 (fa, k 3 )S 12 (fa,k 2 ), (4) 

(3) is the only constraint to the integrability of an open boundary systemEl. Normally, a c-number reflection matrix 
indicates cither a boundary field in the spin chain models or a scalar potential in a fermion system. However, in the 
present model, the boundary impurities have internal degrees of freedom and the spin-exchange processes must be 
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included when an electron is reflected by the boundary. That means the reflection matrix Rj must be an operator 
matrix rather than a c-number one. Below we derive the integrable conditions for the different choices of J and V. 

(i) J = 2, V = 3/2 case. Since the reflection process only consists of one-electron effect, it is convenient to consider 
the single-particle eigenstate. The Schrodinger equation for this case reads 

- $(x + 1) - V(x- 1) = EV(x), Kx<L (5) 
-*(2) = EV(1), x = l (6) 
-V{L-l) + (Vi + J,a -S)^(L) = E^{L), x = L. (7) 

We make the following ansatz for the wave function ty(x) 

V{x)=A + e lkx +A_e- lkx , (8) 

then the eigenvalue E takes the value — 2cosfc. From (6) we obtain 

Similarly, from (7) we obtain 

m-^^-- C + + T + T a m - (10) 

A + e lk + {Vi + Ji(T ■ S) 

The self-consistent condition for A± constrains the value of k by the following eigenvalue problem 

-R{k)A- =e- 2lkL A_. (11) 

For arbitrary AT— particle case, Rj (kj ) must satisfy the reflection equation (3) . It is known that the two-body scattering 
matrix takes the form 

Sd^-m) = - qj - qi+iPjl , (12) 
Qj -Ql-i 

where qj = l/2tan(kj/2), Pji is the spin exchange operator^ 

As shown in appendix A, the reflection matrix Rj(kj) must take the following form 

1 J ' qj -ic + i(S + 1) qj+tc + iS 

where c is an arbitrary constant. By virtue of the relation 

(a 3 • S) 2 = 5(5 + 1) - Uj • S, (14) 
we obtain <jj ■ S has two eigenvalues. Those are 5 and —(5 + 1). For a j ■ S = 5, comparing (10) and (13) we have 

e MU j )_ e^+(V l + SJ t ) 

e-v + iVi + SJiY { ' 

For aj ■ S = —(5 + 1), we have 

q 3 -ic-iS qj +ic- i(S + 1) e ik i + [Vi - (5 + 1) J,] e~ ik i + (V + 5J t ) 



qj + ic + iSqj-ic + iiS+l) e~ ik + [V - (5 + 1) J l ] e ik + (Vi + SJi) 
The solution of the above equation for Vi and Ji leads to the integrable condition 



(16) 



1 „ = 5 2 + 5 + l-(c- l) 2 

(5+i) 2 -C 2 ' 1 (^+l)2_ c 2 



~ rci , 1\9 9> V — ,„ 1,„ n I (1^) 



or 
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s 2 + s + i-{c-\ ) 2 

( S +I)2_(1_ C )2' ( S+ l)2_ (1 _ c) 



^ = 75 i Hi , . 2 ■ ( 18 ) 



Obviously, to ensure the integrability of the present model, Jj and Vi must be related to each other and parameterized 
by a unique constant c as shown in (17) and (18). Notice that (18) is dual to the solution of (17) under the shift 
c — > 1 — c. Without losing generality, we use (17) in the following text. The reflection matrix (10) in the integrable 
case can be rewritten as 



Qj — § Qj — ic + i(a ■ S + 1) qj + ic ■ 
qj + | qj + ic + i(S + 1) qj+ic + iS 



Rj{qj) = ^ ^ . (19) 



For an N— particle system, suppose the wave function initially has an amplitude £ . When the j — th particle moves 
across another particle, it gets an S— matrix Sji{qj — qi). At the right boundary, it is completely reflected back 
and suffered a factor exp(2ikjL)Rj(qj). Then it begins to move toward the left boundary. When it reaches the left 
boundary, it will be kicked back and suffered a factor —1. Finally it arrives at the initial site and finishes a periodic 
motion. Therefore we have the following equation 

— ^j'i-i ' ' ' SjiSji ' ' ' ^-1^+1 ' ' ' SjNfy e2t }L ^jN ' ' ' ^jj+iCo = Co, (20) 

or more neatly 

Sjj-i ' ' ' S ji S ji ' ' ' S ti-i S ?i+i ' ' ' S jN R o S jN ' ■ ■ S jj+i(a = -( 3 _ f ) 2L Co, (21) 

where S% = Sji(qj ± qi). Eq.(21) is just the reflection version of Yang's eigenvalue problerr0. Its solution gives out 
the Bethe ansatz equation. To keep the continuum of the text, we attribute the solution of (21) to appendix B. Below 
we write down the result 

( SL_| )2i+1 = _ q j -i(S + l-c) jj g qj-r^-i g ^■-rA. + i 
Qi + I qj + i(S + l- c) ll J± ?j _ rgj + z ll q . - r \ a -i 



\ a +i{S+±~- c) \ a + i(S - § + c) -j-j- JX- A a - rqj + § TT TT A " - r\p + i 

\ a - i(S + \ - c) A Q - i(5 - | + c) U 1 = 1 A Q - r<& - § " ^ ^ A Q - rA^ - z ' 1 j 

where M < N/2 is the number of spin-down electrons. The eigenenergy is given by 

N 

E = 2N -V — = . (24) 

The derivation of the integrable conditions and the Bethe ansatz equations for other choices of J and V is similar. 
Here we omit the detail and only write down the result. 

(ii) J = 2, V = — 1/2. In this case, the scattering matrix takes the form 

Sil(qi-qi) = Itan- 1 ^. (25) 

qj - qi -i 2 2 



The integrable condition is 

■h , — r- r, = — i; , (2(i! 

and the reflection matrix is 



5' 2 + 5 + l-(c+i) - 
(5 + i) 2 -c 2 ' " l ~ (5 + i) 2 -c 2 



9j + ^ qj — ic — i(a ■ S + 1) qj + ic — ia ■ S 

The Bethe ansatz equations read as 



^.+i c + f(S + i) 'q+ic-iS ' (27) 
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( 1j + j ,2L+i _ Qj -i(S+l + c) -j-r tt Qj - rX a + § 

V qj +i(S+l + c) 11 11 gj - rA a - i ' ^ 

A a - i(5 + \ - c) A a - i(S - | + c) TT TT A « r< ?J § = TT TT A « ~ rA /3 ~ j (2 g) 
X a + l (S+±-c)X a+l (S-± + c) 1 = 1 1J Aq, - rg,- + § 1 = 1 11 A Q — rXp + z ' 

(iii) J = —2, V = 1/2 case. In this case, the two-body scattering matrix takes the form 

^-"'' V-V*? ' *-5*4- (30) 

The integrable condition is the same to that of case (i) and the Bethe ansatz equations are 

( Si + I ^2L+l = Qj + jOg + 1 - c ) TT TT gj ~ rA « ~ § ( o U 

Qj-r ?,-<(5+l-c)ll 11^-rAa + i' 1 J 

A a - i(5 + | - c) A a - i(5 - | + c) -j-j- JX- A a - r gj - - | = „ il A a - rX p - j 

A Q + i(S + \ c) X a + i(S \ + c) 1 = 1 1 = 1 A Q - r 9j + § " 1 = 1 11 A Q — rXp + i ' 1 J 

(iv) J = —2, V = —3/2 case. In this case, the two-body scattering matrix is 

-qi) = , <7j = x tan — . (33) 

Qj — Ql+ l z * 

The integrable condition is the same to that of case (ii) and the Bcthc ansatz equations read as 

( Sj ZJ, X2L+1 = gj - £0g + 1 + c ) TT TT gj -rqi-i y^r gj - rA a + § 

g,+r " %+^+i+c) i = iii gj - rg; +iii g ._ r A a -i' 1 j 

A a + z(5 + | - c) A a + z(S - \ + c) -|-|- JX- A a - rgj + \ = ^ -pt A a - rAg + z 
A Q - i(S + i - c) A a - i(S - i + c) 1 = 1 X a rq 3 f 1 = 1 11 A a - rA^ - i ' 

In the following text, we concentrate ourselves on case (i). Other cases can be treated with the same procedure 
without any difficulty 

III. GROUND STATE PROPERTIES AND LOW-LYING EXCITATIONS 
A. Boundary bound state and zero temperature magnetization 

Unlike in the homogeneous systems, some boundary modes may exist in the present model due to the presence of 
the impurity. From the Bethe ansatz equation (22) we can see an imaginary mode q = i(S + 1 — c) is a possible 
solution as long as c < S + 1. However, such a mode is not always stable in the ground state. As we can read off 
from (24), the real modes form a band —2 < e(q) < 2. The energy of the imaginary mode falls either below the band 
or above it. Therefore, the imaginary q mode defines a true bound state around the impurity. Below we discuss the 
formation of the boundary bound state and zero temperature magnetization for different c values. 

(i)c > S + 1. In this case, the system falls into the regime of ferromagnetic Kondo coupling. There is no boundary 
bound state of charges at low energy scales. For convenience, we put N to be odd in this case. In the ground state, 
all the modes {qj} and {A Q } take real values. Define the quantities 
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M N 

Zt(q) = ^{E E 6 ^ - rA «) "EE ^ - + + 4(g)} + *i(g), (36) 

r— ± a— 1 r— ±j — 1 



Z iW = ^E X>( A ~ r %) -EE ^( A - r\ p ) + 4(A) + 4(A)} 
where O n = 2 arctan(2a;/n) and 



r=±j = l r=±/3=l 



4(g) = 20! (g), 4(g) - 2 (c-s-i)(9), (37) 
05(A) = 2 (A) 0t(A) = 0( 2c+2 s-i)(A) - 0( 2c -2s-i)(A). (38) 



Notice the zero modes are forbidden in an open boundary systemEE Obviously, Z£(gj) = nlj/L, Z 8 L (\ a ) = irJ a /L 
give the logarithmic version of the Bethe ansatz equations, where Ij and J a are integers. In the ground state, both 
{Ij} and { JqJ are closely packed numbers from 1 up to N and M respectively. The cutoffs of q and A are defined as 
ZjXQ) = n(N + 1/2)/ L, Z S L (A) = tt(M + 1/2)/ L. Define the density functions as 

^-h^-k^ (33) 



We have the relation 



1 . r A 1 



pl{q)dq=-N, / p s L (X)dX=-M. (40) 

Q L J-A L 

As demonstrated by many authors0, A — > 00 fqt-the ground state in the thermodynamic limit L — > 00, since any 
hole in the real A— axis induces an excited states. From(40) we deduce that N = 2M, which means the residual 
magnetization is 

M s =S+^N-M = S. (41) 

Therefore, the impurity moment can not be screened any more in this case. 

(ii) 5 + 1/2<c<5'+1. In this case, the system falls still in the ferromagnetic Kondo regime. Due to the stronger 
Kondo coupling and the weaker impurity potential, one boundary bound state occurs with q = i(S + 1 — c) . Since 
I q I < 1/2, we can see the energy carried by this mode 

^^- l-^+l-c)^ - 2 (42) 

is below the band. Therefore this state is much more stable than any other states. Similarly, we can define the 
functions ^£(g) and -Z£(A) for the real modes. In this case, the summation for the real charge modes runs from 1 to 
N — 1 and <j> l c (q) and (j) l s (X) are given by 

9%) = -02(S+l-c) - 02(c-S)(g) - 02(S+2-c)(g), (43) 
4(A) = 0(2c+2S-l)(A) + 0(25+3-2c)(A). (44) 

Notice the phase shifts (jf c s induced by the open boundary are not changed. By integrating p s L (X), we obtain again 
N = 2M. It seems that the local spin is still unscreened. However, we remark the localized electron and the impurity 
form a spin-(S' + 1/2) composite. When c — > 5+1/2 + + , J, V — > — oo, the composite behaves as a perfect local spin 
with spin-(S' + 1/2). It is the composite rather than the original impurity interacting with the host effectively. In 
such a sense, we can say the local moment is partially screenj6d|-Jf we include another half chain interacting with the 
impurity, the problem becomes a two-channel Kondo problerrt3'E3 and we expect the realization of Furusaki-Nagaosa's 
conjecture. 

(iii) — [S — 1/2) < c < S + 1/2. The system falls into the antiferromagnetic Kondo regime. The boundary bound 
state q = i(S + 1 — c) is no longer a stable state since it contains much higher energy. Like in the c > S + 1 case, 
the ground state is described by closely packed real q— modes and A— modes. The residual magnetization is S — 1/2, 
which indicates the Kondo screening, like the case in the traditional Kondo problem. 

(iv) — S < c < — (S — 1/2). Notice that both Ji and Vi are very large in this case. It seems Vi dominates over Ji 
and prevents the conduction electrons near to the impurity. The residual magnetization is S, which means there is 
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no Kondo screening. Such a situation is very different from that of a Kondo problem in a 3D Fermi liquid, where 
the impurity potential does not change the fixed point and only induces the renormalization of the Kondo coupling 
constantEa. Such a phenomenon strongly suggests the charge-spin cooperation plays an important role in the ID 
Kondo problem. 

(y)—(S + 1/2) < c < —S. In this case, the system is still in the regime of antiferromagnetic Kondo coupling. From 
the Bcthc ansatz equations(22-23) we can see an imaginary spin mode A = i(S — 1/2 + c) may assist the formation 
of a boundary bound state with q = i(S + c) . This mode carries much lower energy than that of a real mode and is 
therefore stable in the ground state. The functions Z^q) and Z S L (\) can be defined in a similar way. Notice here the 
summation for the real charge modes runs from 1 to N — 1 and that for the spin modes runs from 1 to M — 1. 4> l c (q) 
and <pl (A) are defined as 

4>i(q) = -02(x+s-c)(q) - 2{s +c+i)(q) - 2 \ S + C \{q), (45) 

<AW = 0(2S-2c+l) - 9\2S+2c-3\W- (46) 

In the thermodynamic limit, we have N = 2M — 1 in the ground state. Therefore, the residual magnetization takes 
the value of S — 1/2, which indicates a typical Kondo screening. In fact, both the positive Jj and negative Vi assist 
one conduction electron to form tight-bonding pair with the local moment. The bounded spin mode is a signal that 
the local spin-(5 — 1/2) pair is formed. 

(vi)c < — (S + 1/2). In this case, both the ferromagnetic Kondo coupling (Jj > 0) and the repulsive impurity 
potential forbid the electrons to form a bound state with the impurity. Z^(q) and Z S L (X) can be defined in the same 
forms as (36) but with 

#(s) = -0a|sr +1 _ c |(g), (47) 
#W = 0(2S-2c+l)(A) - 0|2S+2c-l|(A). (48) 
The residual magnetization is still S and no Kondo screening occurs. 

B. Low- lying excitations 

As we discussed above, there is a finite gap between the local level of the boundary bound state and the conduction 
band. At low temperatures, the excitation of an electron from the local level to the Fermi surface must overcome 
the energy gap and is thus unimportant to the low temperature thermodynamics. In this subsection,.-we discuss the 
low-lying excitations near the Fepni surface. Notice that the current excitations in a periodic systemoa are forbidden 
here due to the open boundaries^! 

In the thermodynamic limit, we denote the density functions as p c (q) , p s (X). The energy per unit length is thus 



e 



e° c (q)p c (q)dq, (49) 
-Q 



where e®(q) — 2 — [q 2 + (1/2) 2 ] 1 and p c (q) satisfies the following integral equations 



p c ( q ) = a 1 (q) + ^ z [^-cf>' c (q)-6(q)}- I a 2 i q - f /),vU/ )>h/ + / "i (</ - A)/i.,( Ai</A. (50) 



A 
-A 



PsW - ^[^(A)-<5(A)]+ J a 1 (X-q) Pc (q)dq- J ^a 2 (X- X')p s (X')d\' . (51) 
Here axe the derivatives of (f>^ s + 4> l c s ; a n (q) = 9' n {q)/2i:. By substituting (50) into (49), we readily obtain that 

Q e c (g)oi( g )dg + -U C e c (q)cj>' c (q)dq + f e s (X)^(X)dX - 2^°(0)], (52) 
Q 4ttL J_ q J_ a 

where the first term is just the energy of the bulk and the second term is the energy induced by the impurity and the 
boundaries; 

e c (q) = e°M - [ Q a 2 (q - q')e c (q')dq' + a x (q - X)e s (X)dX, (53) 

J-Q J-A 
rQ /-A 
e s (X)= ai (X-q)e c (q)dq~ a 2 (X - X')e s (X')dX' (54) 

J-Q J-A 
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are the dressed energy of the systemS0. It seems that the impurity and the open boundary have no effect to the 
dressed energies up to the order 0(1/ L). If the system is at an excited state, some particles in the lower Ij- and 
J Q -states will be moved to higher Ij- and J Q -states and leave some holes in the Fermi seas. We denote the momenta 
of excited particles as q p and A p , while those of holes as qh and A^. The bare change of the density distributions p c (q) 
and p s (A) are thus 



r=± p h 
r=± p h 

In addition, <5p^ and 5p b s have feedback effect to p c and ,o s in the Fermi seas and induce the change of p c and p s 
Spc{q) = EE a ^(l- r 1h) - ^2a 2 (q-rq p ) + E a i(<? " r\ p ) - E«i( A ~ rA/,)] 

r=± ft p p h 

a 2 {q - q')Sp c (q')dq' + / a x (q - A)<5p s (A)dA, (57) 

<5^(a) = i e E ai ( A ~ r? p) ~ E ai ( A ~ r " 9fe ) ~ E a2 ^ A ~ rA p) + E a ~ 2 ( A ~ rA/i )] 

r=± p /i p /i 

+ / a 1 (\-q)Sp c (q)dq- a 2 (\ - \')5p s (\')d\' . (58) 

J-Q J-A 

The excitation energy is 

AE = 4(l P ) - E + [ Q e°MSPc(q) 

p ft 

= E e ^ - E £c (*) + E £s ( a p) - E e M- ( 59 ) 

p h p h 

It is just the same expression for those of a pure systemHS. That means the impurity does not affect the bulk 
excitations up to the order 0(1/ L). The. dressed energies e CjS thus can be treated as quasi-particle energies in the 
so-called Landau-Luttinger liquid thcoryiEI. Nevertheless, the impurity indeed affects the thermodynamics in the order 
of 0(1/L) since it induces the changes of the_density of states in this order as shown in (50), (51). We remark there 
are also string excitations in the spin sectorEl The excitation energy can be-expressed as the same form of (59) and 
the strings induce only different scattering matrices from those of spin holes 
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IV. THERMODYNAMICS 

In this section, we derive the thermodynamic equations of the present model via thermal Bethe ansatzEG. We 
shall omit the excitations which break the boundary bound state since the excitations are accompanied by finite 
energy gaps and their contributions to the low temperature thermodynamic quantities are exponentially small. 

A. Thermodynamic Bethe ansatz 

At finite temperatures, the solution of the Bethe ansatz equations are described by a sequence of real {qj} and a 
variety of {A a } strings. From the Bethe ansatz equations we obtain 



11 °° 

M) + p h c(q) = + 4^#(?) - 2l%) " MpM + E (60) 

rn—1 

oo 

Pa,m(A) = j-r^mW - ^(A) + [m]p c (A) - E A mn p s ,nW, (61) 
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where p Sy7n (p% m ) are the densities of the m-string (holes) ;[n] is an integral operator with the kernel a n (q); A mn = 
[m + n] + 2[m + n — 2] + ■ • • + [\m — n\ + 1]; <j) s ,m — Sjli fisi^ + i((jn + l)/2 — j)]. In a magnetic field H, the free 
energy of the system can be written as 

F/L = f [e° c (q) - n - ^] Pc {q)dq + J2 nH [ Ps,nWdX 
J n=l •* 

-T J [{p c + p h c ) ln(p c + p h c ) - p c \np c - p h c \np h c ]dq 

CO „ 

~ T E / ttP*> n + Pin) HPs.n + Ps, n ) - Ps,n In Ps,n ~ Ps, n ^ Ps,n\dX. (62) 



n=l ' 



At equilibrium state, the free energy must be minimum. This gives 

In = £ " ~ M T ~ ^ + [2] ln(l + r,- 1 ) - ^[n] ln(l + C" 1 ), (63) 

n 

ln(l + C„) = ^ - [n] ln(l + v^) + E A ™« Ml + &)> (64) 

m 

where 77(g) = p^(q)/ p c (q), £ n (A) = p% n (X)/p Sjn (X). The above equations can be rewritten as 

lnr, = + {[2] - [l]G}ln(l + r?" 1 ) - Gln(l + Ci), (65) 

lnC„ = G[ln(l + C„-i)+ln(l + Cn+i)], n > 1, (66) 
InCi = -Gln(l + r y - 1 ) + Gln(l + C 2 ), (67) 

where G is an integral operator with the kernel [2 cosh(7rA)] _1 . The free energy takes the form 

F = Fa + —Fimp + j^F e dg, (68) 



Fo/i = -TJ ai(q) ln(l + ly" 1 (<*)]<*?> (69) 
F edg = ~Tln(l +t?- 1 (0)] + ^T^ln[l + C^O)] 

n 

/ #'(g) ln(l + rr^dq ~ f E / <^«'( A ) ln ( X + C')^, (70) 

Fimp = ~J- f € '(?) + T 1 )^ - £ E / <^«'( A ) ln ^ + C 1 ) A ( 71 ) 

where Fo is the free energy of the bulk, F e d s is the free energy of the open boundary and Fi mp is the free energy of 
the impurity which contains a part of the boundary energytj, the charge-fluctuation term and the contribution of the 
impurity spin. 

B. Residual Entropy 

To give a detailed picture of the ground state configuration, we study the residual entropy of the impurity. When 
T — ► 0, rj — > exp[(e c — /i)/T], for & . < a. Therefore the driving term in (67) tends to — oo. That means £i — > and 
all the other ( n take constant valueaHaS with 

+ sinh 2 (na;o) 1 H 

Ct = — T2 1 > x o = (72) 

sinn xo ^ J 

Here we consider the spin part of the impurity free energy F? mp , which is the only relevant term to the residual 
entropy and the magnetization. With different c, the ground state configurations are very different. Below we derive 
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the residual entropy for different parameter regions. Since at zero temperature, £n are variable- independent, the 
integral operators [n] are equivalent to unit and a n arc equivalent to (5(A). This makes the calculations very simple. 
For half-integer or integer c, the impurity spin is split into two effective spins Si and S2 and the residual entropy 
reads 

Sres = \sgniS,) ln(2|^|) + \sgn{S 2 ) ln(2|S 2 |). (73) 

The following discussion is for non-half-integer c case. 
(i)c > S + 1. In this case, 



1 25 

2^8, n W —* a n,2 Cj +2S-l(A) - a„,2 Cj -2S-l (A) - a Q (A) ^(5„,2c / -2S+2;-2, (74) 



1=1 



where c/ is defined as: if n/2 < |c| < (n + l)/2, cj = \n\/2; a = 2(\c\ — cj); and 



minijn.n) 

^ a n +l+m-2k- (75) 



fe=l 

For T -» 0, 

1 + C+ 1 25 



- 4 Tln i + ^ +M " 1 + ^ T E ln (! + C^-25 +2; -2 _1 )- (76) 



+ C 2cj -2S-1 " ;=1 

Therefore, the residual entropy takes the value 



1 2c/ + 25-l ^, 2(cj-S-l + /) 



S r p.s — ~ In 



i-Vln - ^---^^> . (77) 



2 2cj-25-l ^ v /4(c7 - S - 1 + I) 2 - 1 

The magnetization at zero temperature is M = S, which is the same result given in sec. III. 

(ii)S + 1 > c > S + 1/2. In this case, a boundary bound state exists in the ground state and cj =5 + 1/2. 



2tt 

The residual entropy reads 



I 2S-1 
-^-4>l, n W ~> fl ",4S + ffln,2 ~ a a 5 "-,2J+l - (78) 



2,5 1 2Z + 1 

S res = ln(2 V2S) ln ' < 79 ) 



(iii) 1/2 < c < S + 1/2. No bound state exists and 

2c/-l 

.A* 'CAI 



9^^s,n W ~~ * a n,2S+2cj-l + «<n,2S-2c 7 + l ~ Aq ^n,2S-2 C/ +2/ • (80) 

This gives 



2cj — 1 

(iv)c/ = 0. In this case, c has no effect to the residual entropy and 

^<„'(A) - 2a„, 2S . (82) 

The residual entropy reads 
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S res = ln(2S), (83) 

a similar result to that in the conventional Kondo problem@. 
(v)-(S- 1/2) < c < -1/2. 

0a,n M — * a ",2S+l+2 C/ + On,25-X-2c/ — O q ^ £n,2S-2 Cj +2i, (84) 



and 

2cj 



(vi)-5 < c < -(5" - 1/2). In this case, a = S - 1/2 and 

2S-1 

rA 4 r A^l — > «.„ ,1 c — n.„ 

2tt 



j 2S-1 

I - 0s n W — * a ",4S - a a <5«.2Z+1- (86) 

Caution should be taken in this case since ln(l + Ci ) appears in the free energy, which is divergent when T — ► 0. 
Note that 21n(i - ln(l + r/ _1 ) + ln(l + ( 2 ) as shown in (67) and ln(l + ?7 _1 ) takes the asymptotic form AT^ 1 + BT. 
Therefore rj contributes nothing to the entropy. The residual entropy in this case reads 

S res = ln(2 V2S) " E ln ■ ( 8? ) 

(vii)— (S + 1/2) < c < —S. In this case, boundary bound states exist in the ground state and cj = S. 

2S-1 

j) 1 ( X) — KL4C11 — n„ i — a,. 

2tt 



^ 2S-1 
7T~ ti W — * a n,4S+l — dn,3 — da > , <5 n .2+2i- 

/7T ' * • 



grM = lln^±l- Vln - 2 ^ + 1 ) . (89) 
2 3 f- ,/(4(Z + l) 2 -1 V ' 



, , vW+T 

(viii)c < — (5 + 1/2). The situation of this case is very similar to that of case (i) and 

2S 

th I A I — » flu o c j o*> _ i i — n„ o o i i — n~. 

2tt 



^ 2S 

^^s.n W ~^ an,2S+2c, + l - a n ,2 CI -2S+l — Q ^ ^n,2cj-2S+2i- (90) 

1=1 

The residual entropy reads 

2S 

s r „ - h VP. + - - g m V4 g:-/ + y _ , . 

From the above discussion we can see that the impurity-bulk coupling not only changes the effective energy scale as 
in the conventional Kondo problem, but also renormalize the effective strength of the impurity spin. Generally, the 
coupling splits the impurity spin into two effective spins (ghost spins). At low temperatures, the local spin configuration 
is strongly interaction-dependent rather than simply screened or decoupled as understood in the conventional Kondo 
problem. This strongly suggests that the correlation in the bulk has a non-trivial effect to the low-temperature 
behavior of the impurity. 
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C. Local Landau-Luttinger liquid description 



Despite the finite residual entropy in the ground state, the leading order of the low temperature specific heat 
should be T-linear since there is no extra degrees of freedom to induce the overscreening effect. Therefore, Nozieres' 
Local Fermi liquid theoryLJ can be used with a slight modification. Consider a non-interacting ID open chain. The 
"momenta" k' is quantized as 

k' = k+^-6 Q (k), (92) 

where k — irn/L (n positive integer) are the "momenta" of the pure open boundary system and So(k) is the phase 
shift due to the particle-impurity scattering. The change of the density of states is therefore 

{D <*> - hWi- (93) 

where e(k) is the single-particle energy. Fop-the interacting ID electron systems, the quasi-particles can be defined 
in the charge- and spin-sectors respectively^. The phase shifts are generally functionals of the distribution of the 
quasi-particles 5 c [n c (q),n s (X)], S s [n c (q),n a (X)] and can be expressed as 

Sc(q) = S c>0 (q) + Qcc{q,rq')5n c {q') + £ £ 6 ca {q, rX)Sn s (X), (94) 

r=± q'=j£q r=± A 

S S (X) = 6 S>0 (\) + ^ ^sc(A, rq)Sn c (q) + £ ]T 6 ss (X,rX')6n s {X'), (95) 

r=± q r=±A'^A 

where (5 Cj o and S s q are the bare phase shifts induced by the impurity; 0's are the phase shifts of the particle-particle 
scattering; 5n c , s — n c,s ~ n c s is the change of the quasi-particle distributions induced by the impurity. In another 
hand, Sn CiS in the ground state take the form 

Sn c (q) = Sn c (X) = (96) 

The equations (94) and (95) are thus reduced in the thermodynamic limit L — > oo to 

i rqp i r^F 

S'M = <o(3) + 7T / 6> cc (q,q'KWW + 7r 6' cs (q, X)S'(X)dX, (97) 

27r J-q F 27T J-\ F 

i rlF i r^F 

S' s (X) = 5' sfi (X) + — e sc (X,q)5' c (q)dq + Tr C (A, X')8' s (X')dX', (98) 

where qp and Xf are the Fermi "momenta". Notice o(l/L) terms have been omitted in the above equations. The 
low-temperature thermodynamics of the impurity is therefore charactered by two constants 

5D c {q F ) = - — 77 — r, dD s {X F ) = - — 77--^, 99) 
2ne' c (q F ) 2ir e' s {X F ) 

where e CjS are the quasi-particle energies. 
In our case, 

*co(g) = Ml), *.,o(A) = ^s(A), (100) 
9' cc (q,q') = -2wa 2 (q-q'), 0' cs (q, A) = 2wai(q — A), (101) 
9' sc (X,q)=27ra 1 (X-q), 8' SS (X - A') = -2na 2 (X - A'). (102) 

Therefore, 

S'M = <f>'M - [ Q *2(q - q')S' c (q'W + [ ai(q - X)6' s (X)dX, (103) 
J-Q J— A 

rQ /-A 

5 S (A)=<^(A)+ / 01 (A - q)8' c (q)dq a 2 (X-X')6' s (X')dX'. (104) 

J-Q J-A 
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The variations of the densities of states at the Fermi surfaces induced by the impurity are given by (99) with e CiS the 
dressed energies (53-54). For the low-temperature specific heat, the following relation holds 

SC _ 8' C (Q) 8' S (A) 

Co Pc(Q) Ps(AY [ ' 

where SC is the specific heat induced by the impurity and the open boundaries and Co is the specific heat of the bulk 
(per unit length); p c (q) and /O s (A) are the quasi-particle distributions in the ground state. 

V. CONCLUDING REMARKS 

In conclusion, we propose an integrable model of a boundary impurity spin coupled with the integrable open t — J 
chains. In our model, The "fine tuned" effect in the periodic integrable models is overcome and the interaction term 
takes a very simple form. The coupling constant J can take arbitrary value without destroying the intcgrability of the 
Hamiltonian, while in the periodic models, there is a constraint to J. Though a similar impurity can be introduced 
in the periodic models, the parameter c must be real (imaginary in our case) which describes a weak-linked impurity 
to the bulk. The interaction only affects the energy scale (Kondo temperature) but does not change the fixed point 
of the system. With an imaginary c, the model Hamiltonians constructed for bulk impurities are non-Hcrmitian and 
their spectra generally lie in the complex plane rather than in the real axis completely. In our model, real c defines 
a Hcrmitian Hamiltonian due to the reflection symmetry, and the coupling constant J meets all physical situations. 
Some new phenomena driven by the impurity-bulk coupling have been found, which can never appear in the periodic 
models as well as in the conventional Kondo problems: (i)The stronger coupling J may split the impurity spin into 
effective "ghost spins" S—cj + 1/2 and S+Ci — 1/2. That means the coupling not only changes the energy scales (Kondo 
temperature) as in the conventional Kondo problem but also renormalizes the effective strength of the impurity spin. 
It seems that the strength of these ghost spins does not change via temperature. Such a phenomenon reveals a pure 
correlation effect. (ii)Depcnding on the strength of the coupling, the system may show behavior differing from those 
of the conventional Kondo problems. A typical example is that the scalar potential may destruct the Kondo screening 
even in the antiferromagnetic regime. (iii)Thc residual entropy is strongly coupling-dependent which hints the local 
spin configuration near the impurity is very complicated rather than simply screened or decoupled as understood in 
the conventional Kondo problem. 

One of the authors (YW) acknowledges the financial aids of AvH-Stiftung and CNFNS. 

Appendix A 

Since the model is spin SU (2) invariant, the reflection matrix must take the following form 

Rj{Qi) = f{QiMQi) + i^-S], (A.l) 
where f(q) is an arbitrary function and g(q) is determined by the reflection equation. With an S— matrix 

Sjl( qj ,qi)- q ._ qi+h . (^-2) 

from the reflection equation (3) we have 

[<u - m + hPji] [gin) + iaj ■ s] [qj + qi + hPji] [g{qi) + m • s] 

= [g(qi) + im ■ S][ qj +qi+ hPjlMqj) + itTj ■ S][qj - qi + hPji}. (A3) 

By virtue of the relation (a ■ S) 2 = S(S + 1) — a ■ S, we obtain 

j qj [2g{ qj ) -i}- q] = jqi[2g(qi) - i] - qf. (A4) 

Since qj,qi are arbitrary, the left hand side and the right hand side of the above equation must be a constant C. 
Therefore, 

, \ q c i , . 

g(q) = - + — + -. A5 
7 79 2 
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Substituting (A. 5) into (A.l), we readily obtain 

R { ) = eMa) g, ~ + s +1) £ + ^ ± s , (a.6) 

JV JJ qj - icy + ij(S + 1) + v 7 

where <£>((/) is an arbitrary function and c is an arbitrary constant. 

Appendix B 

In dealing with the open boundary integrable models, we often encounter the following eigenvalue problem 

S n-i ■ ■ ■ S Ii S ti ■ ■ ■ S ±i-l S ^+l ■ ■ ■ S tN R i S lN ■ ■ ■ Sj~ j+ ii>o = X^ = e fo)V>o, (3.1) 

where and Rj are the scattering matrix and the reflection matrix respectively. In our case, up to constants they 
take the forms 

si - gj± + g ' + : 7 .^ (b.2) 

J qj ±qi+vy 



qj — icy + iy(<Jj ■ S + 1) qj + icy + iy<Jj ■ S 
qj — icy + iy(S + 1) qj + icy + iyS 



^ _ V j - / -r 6 n vj / ~r~ / ~r~ 1 j {B 3) 



Define V->o = {Sjj-i ■ 1 1 Sji) 1- 0o- can be rewritten as 

Sjf' Sjj-iS jj+ i ■ ■ ■ S+ N RjSj N ■ ■ ■ Sjj^Sjj^ ■ ■ ■ S^ipa = X^o = e(qj)^ - (BA) 
For convenience, we introduce an auxiliary space r and define 

U T (q) = S T jS Tl ■ ■ ■ S T j_! ■ ■ ■ S^ N R T aS TN ■ ■ ■ S T j +1 S r j_ 1 ■ ■ ■ S t1 S t j, (B-5) 
with q T = q. Obviously, S~^(qj) = P T j and 



Since satisfy the Yang-Baxter relation 



S- T ,(q - q')S±{q ± qj)S±j(q' ± Qj ) = S^q' ± qj )S±(q ± ^^(q - q'), (B.7) 
from (3) we can easily show that U T (q) satisfies the reflection equation 

S; T ,(q - q')U T (q)S+ T ,(q + q')U T ,{q') = U T ,(q')S+ T ,(q + q')U T (q)S^ T ,(q - q'). (B.8) 
Therefore, the eigenvalue problem (B.3) is reduced to Sklyanin's eigenvalue problemEl. Furthermore, we define 

T T (q) = S To (q)S~ N (q) ■ ■ ■ S~ j+1 ■ ■ ■ S'^ ■ ■ ■ S^S^(q), (B.9) 

where 

Sro = g+ : 7c+ : 7 % s - (b.io) 

q + lyc + lyb 

Then 

U T (q)=T- 1 ( -q)T T (q), (5.11) 

and T T (q) satisfies 

S- T ,(q - q')T T {q)T T ,{q>) = T T ,(q')T T (q)S; r ,(q - q'). (5.12) 
For convenience, we introduce the notations 
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Notice that 



Therefore 



UAq) ~{c(q), D(q))> 1Aq) -{c(q), d(q) ) ' 



T2S-/(-q i!)r 2 = S+(q)(l + — >-), 

q + qj 

r 2 5 T0 {-q z 7 )r 2 = S T0 (q)[ q _ _ + (5.13) 



rp-if \ q-ici- ij{S +1) tt q + qj _, . . ,„ > 

a-icy + iMS + 1) J -- L o + qh + i-i 



q-icj + i^{S + 1) q + qj + i 7 
From (B.7) and (B.12) we have the following commutation relations 

[c(g), %')] = WW?) - d(q)a(q')}, (B.15) 

[A(q)+D(q),A(q')+D(q')} = 0, (£.16) 
A(q)B( P ) = {q + P \ { r P ~ ll \ Btv)A{q) - %- -B( P )D( P ) 

+ 7 x-B(pM(«z), (B.17) 

5(g)B(p) = b- p + W« + > + W B(p)D{ q ) - . 2 ^ + y l^Dfr) 
q - p){q + p + n) (q-p)(2p + n) 

+ - l!2%±^ B( gM(p) , (£.18) 
(2p + z 7 )(g + p + i 7 ) w VF; ' v ; 

where 5(g) = (2g + i^)D{q) — i^A(q). To construct the eigenstates of Xj, we define the pseudovacuum |fi > as all 
spin up. From the definition we know 

a(q)\n>=\fl>, c(q)\n>=0, 
d(q)\n >= l±mizill JT _J^_| fi > . (5.19) 

Therefore 

A(g)|Q>=|Q>, C(g)|fi>=0. (B.20) 

Notice that 
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With the commutation relations (B.15) and (B.19), we obtain 



S(«)|n >= 2 q q - ll( : S Q + c -]\ q -^ s - c \ n i {q + q ^ q - q > } . j n > . (5.23) 

The operator B(q) can be treated as the creation operator of the eigenstate. An eigenstate of X'(q) with M spin 
flipped can be written as n^Li B(A a )\Q >= Q M \^ >• From (B.17) and (£.18) we have the following relations 



M 



A(q)Q M = Q M A(q) Y[ f± 



(q + A a )(q- A a -ij) 



=i (q - K){q + A a + i-t) 



M 



V l l /V a i IT ( A o-A/3 + ^)(Aq+A /3 + 2z7) 

^ i (g + A Q + J 7)(2A ct + i7 )^ M 1 q; 11 (A a -A /3 )(A a + A /9 + i 7 ) 

y 2 ^A a a . , ■ tt (A a - Aff)(A a + A,3 - 27) 

+ («-A a )(2A a + i 7 ) g ^ lAa) 11 (A Q -A,)(A Q + A, + z 7 )' ( } 

d( 3 )Qm - qm^w n {q _ Aa){q+Aa+n) 



_ ^ 2i 7 ( g + ij) - j-r (A a -A^ + z 7 )(A a +A ;3 + 2i7) 

^ (« - A Q )(2A Q + yy) ^M^a) U (Aq _ Wq + Ap + - 7) 

, ^ 4i7A a (g + ^7) (A a - A p )(A a + Ag - ij) 

+ (q + A a + ry) (2A Q + *y) y ^ (Aa) 11 (A a - A )(A a + A, + i 7 ) ' 

where Q%j = B{A\) ■ ■ ■ B(A a _i)B(q)B(A a+ i) ■ ■ ■ B(Am)- Obviously, when we apply A(q) and D{q) to the state 
Qm\& >, two types of terms appear. One is the wanted terms with Qm, the other is unwanted terms with Q%- With 
the relation Xj(qj) = [2(qj + i'f)]~ 1 D(qj) + A(qj), we obtain 

<v) = n gj+ ! agj "! a '^ - (^-26) 

q=i% - A a <&+A a +t7 
The cancellation of the unwanted terms gives the constriction of the parameters A a 

A a - i~/(S + c - 1) A a - i 7 (S - c) (A a + 9j )(A a - % ) 

'-c+ll 11 



A a + i~f(S + c) A Q + ij(S - c + 1) 11 (A Q + ^ + i-y){A a - q 3 + ij) 



{A a + Ap)(A a - Ap - i-y) 
Q +A /3 + 2i7)(A Q -A /3 + 

By replacing A a with X a — i-f/2, we readily obtain the Bethe ansatz equations. 



(A a + A /3 + 2t7)(A Q -A /J + i7)' 
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